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1.0  INTRODUCTION 


The  ability  of  the  FAA  Beacon  Collision  Avoidance  System  (BCAS)  to  detect 
and  properly  respond  to  collision  hazards  depends  upon  the  accurate  estimation 
of  aircraft  vertical  rates.  These  estimates  are  based  upon  altitude  reports 
derived  from  encoding  barometric  altimeters  and  are  quantized  in  100-foot 
altitude  increments. 

In  early  BCAS  testing  it  was  noted  that  the  response  of  the  altitude 
tracking  algorithm  to  a  single  isolated  transition  from  one  quantization  level 
to  the  next  was  a  rate  estimate  of  substantial  magnitude.  Since  Isolated 
transitions  often  occur  when  the  actual  rate  is  negligible,  the  resulting 
error  in  rate  estimation  could  lead  to  Improper  selection  of  avoidance 
maneuver  directions. 

Fus  ther  study  revealed  that  a  large  part  of  the  problem  was  due  to  the 
manner  in  which  the  linear  recursive  "alpha-beta”  tracking  algorithm  responded 
to  quantized  data  inputs..  Reductions  in  tracking  gain*  reduced  tracker 
response  to  isolated  altitude  transitions.  Although  such  gain  reductions 
eliminated  the  over-response  problem,  they  also  substantially  reduced  the 
capability  of  the  tracker  to  respond  to  actual  rate  changes.  As  a 
consequence,  Lincoln  Laboratory  undertook  an  investigation  of  alternative 
tracking  techniques  which  suppress  response  to  isolated  altitude  quantization 
transitions  while  responding  promptly  to  actual  rate  changes.  This  effort 
yielded  improved  tracking  algorithms  for  use  in  collision  avoidance  systems. 
It  also  provided  theoretical  Insight  into  the  general  problem  of  rate  tracking 
with  coarsely  quantized  inputs. 

As  employed  in  this  document  the  term  ’’coarsely  quantized”  refers  to  a 
system  in  which  the  rates  of  interest  are  such  that  at  most  one  quantization 
level  is  crossed  between  samples.  For  the  100-foot  quantization  levels  and 
1-Hertz  update  rate  of  BCAS,  all  rates  below  the  nominal  design  limit  of  6000 
FPM  result  in  coarsely  quantized  rate  tracking  behavior.  If  altitude  tracking 
is  based  upon  ground-based  sensor  data,  the  update  interval  may  be  4  seconds 
or  greater.  Higher  altitude  rates  may  then  result  in  several  quantization 
levels  being  crossed  between  samples  and  the  ability  of  the  tracker  to 
estimate  rate  is  limited  more  by  the  sampling  rate  than  by  the  altitude 
quantization.  With  4-second  update  intervals  it  is  more  appropriate  to  refer 
to  the  tracking  process  as  "coarsely  sampled”.  The  design  principles  used  in 
the  coarsely  quantized  tracker  apply  with  only  slight  modification  to  the 
coarsely  sampled  regime,  as  will  be  described  in  Section  8.0. 


•Broste,  Nels  A.,  "A  Vertical  Tracker  Redesign  for  Active  BCAS,” 
MTR-79W00431,  The  MITRE  Corporation  (March  1980). 
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2.0  MODEL  OF  THE  ALTITUDE  TRACKING  PROCESS 


The  design  of  estimation  algorithms  requires  that  a  mathematical  model  of 
the  measurement  process  be  defined  in  order  to  provide  the  basis  for  both 
mathematical  analysis  and  computer  simulation  efforts.  This  section  discusses 
those  characteristics  of  air  traffic  end  altimetry  which  are  relevant  to  the 
selection  of  an  appropriate  model 

2.1  Error  Characteristics  of  Altitude  Reports 

Observation  of  barometric  altimetry  data  indicates  that  errors  in 
pressure  transducer  output  do  not  vary  significantly  from  sample  to  sample. 
Such  constant  errors  can  be  viewed  as  biases  which  contribute  errors  In 
position  determination  but  do  not  influence  the  accuracy  of  altitude  rate 
estimation.  Hence  minimization  of  error  due  to  sample-to-sample  jitter  is  not 
a  primary  problem  in  altitude  tracking  (although  tracking  should  be  designed 
to  smooth  such  errors  whenever  they  occur).  The  principal  challenge  of 
altitude  tracking  is  to  produce  accurate  estimates  in  the  presence  of 
measurement  quantization  and  finite  sampling  rate. 

2.2  Relative  Significance  of  Rate  and  Position  Errors 

In  collision  avoidance  applications,  the  quantity  of  interest  is  the 
aircraft  position  at  some  future  time  (such  as  the  time  of  closest  approach  or 
the  time  at  which  a  response  to  collision  avoidance  Instructions  would  begin). 
Altitude  projections  are  obtained  by  a  simple  linear  projection  of  aircraft 
motion  according  to 


l 

t 


2  .  £  t 

0  +  0 


(2.1) 


where  zt  is  the  projected  altitude  at  time  t  and  zq  ,  zq  are  the  current 
estimated  position  and  rate.  The  rate  is  assumed  to  be  constant  for  the 
duration  of  the  projection.  If  the  errors  in  these  quantities  are  denoted  by 
,  eg  and  ej  then 
tO  0 

e2  "  e2  +  e*  t  (2.2) 

too 


It  can  be  seen  that  the  longer  the  projection  time,  the  greater  is  the 
significance  of  the  rate  error  in  comparision  to  the  position  error.  For  30 
seconds  projection  time,  a  rate  error  of  200  FPM  is  equivalent  (in  terms  of 
projection  error)  to  a  position  error  of  100  feet. 
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Quantized  altitude  measurements  include  a  quantization  error  which  lies 
between  -q/2  and  q/2.  Any  tracking  which  goes  beyond  mere  uae  of  the  raw 
report  involves  attempting  to  determine  exactly  where  within  the  quantization 
interval  (of  width  q)  the  actual  position  lies.  Hence  for  q  -  100  feet  the 
maximum  reduction  in  projection  error  which  can  be  achieved  by  improved 
position  estimation  is  less  than  or  equal  to  50  feet.  It  can  now  be  seen  that 
position  estimation  errors  are  unlikely  to  be  a  critical  issue  in  the 
selection  of  an  altitude  tracking  technique.  This  is  true  first  because 
errors  on  the  order  of  50  feet  are  not  large  enough  to  significantly  affect 
the  performance  of  collision  avoidance  systems.  Secondly,  experience 
indicates  that  the  tracking  design  cannot  be  expected  to  significantly  reduce 
altitude  determination  errors  below  ±q/2. 

Rate  estimation  errors  however  can  be  critical  to  resolution  success  and 
can  vary  significantly  depending  upon  tracking  technique.  An  error  of  1000 
FPM  projected  over  a  30  second  interval  produces  an  error  in  projected 
altitude  of  500  feet.  Such  errors  are  comparable  in  magnitude  to  the  amount 
of  altitude  deviation  which  can  be  effected  by  resolution  commands.  Hence, 
rate  errors  can  and  do  influence  the  success  of  resolution. 


2.3  Characteristics  of  Aircraft  Trajectories 

Several  aspects  of  tracking  algorithm  design  are  influenced  by  the  types 
of  altitude  trajectories  which  are  to  be  expected  in  operation.  The  following 
paragraphs  discuss  some  of  the  significant  characteristics  of  altitude 
trajectories. 

Pilots  generally  attempt  to  maintain  a  near-zero  altitude  rate  at  the 
desired  flight  altitude  or  to  hold  a  constant  non-zero  rate  in  transitioning 
between  altitudes.  In  some  cases  larger  altitude  changes  occur  in  a  number  of 
steps  as  air  traffic  control  issues  clearances  to  successive  altitude  levels. 
Typical  aircraft  trajectories  might  be  envisioned  as  a  series  of  constant  rate 
segments  with  periods  of  acceleration  occurring  whenever  the  rate  changes.  In 
reality  of  course,  aircraft  never  maintain  exactly  constant  rates.  And  even 
when  attempting  to  hold  constant  altitude  there  is  some  oscillation  about  the 
desired  altitude.  Hence  isolated  transitions  between  adjacent  quantization 
levels  may  be  observed  even  for  aircraft  in  nominally  level  flight. 

Lower  performance  aircraft  tend  to  climb  and  descend  at  rates  between 
500  FPM  and  1000  FPM.  Higher  performance  aircraft,  such  as  jet  transports, 
use  higher  rates,  but  seldom  exceed  6000  FPM.  The  BCAS  system  is  designed  to 
provide  nominal  performance  for  rates  up  to  6000  FPM. 


In  going  from  one  altitude  rate  to  another,  aircraft  are  usually  assumed 
to  accelerate  at  values  between  0.10g  and  0.25g.  This  level  of  acceleration 
allows  the  rate  to  change  substantially  between  quantization  levels.  Consider 
for  instance  an  aircraft  which  accelerates  with  constant  acceleration  from  an 
Initial  vertical  rate  £q  to  a  rate  if.  The  altitude  is  given  by: 


z(t) 


z0  +  z0  t  + 


0  <  t  < 


a 


(£f  -  £q)2  #  4f  ~  40 


z0  ~  -  +  zf  t  t  > - 

2a  a 


(2.3) 


where  a  is  the  acceleration  and  Zq  is  the  initial  altitude.  At  e  clme  the 
final  rate  is  achieved  the  aircraft  position  is 

2  2 
if  -  zQ 

Zq  +  -  (2.4) 

2a 

It  can  be  seen  from  the  above  expression  that  an  aircraft  which  initiates 
a  0.25g  vertical  acceleration  from  level  flight  can  achieve  a  rate  of  more 
than  2400  FPM  before  moving  100  feet  vertically.  This  indicates  that  the 
coarseness  of  the  altitude  quantization  is  insufficient  for  accurate  rate 
tracking  during  periods  of  acceleration.  However  It  should  also  be  noted  that 
periods  of  acceleration  do  not  persist  for  more  than  a  few  seconds.  At  0.25g 
for  instance,  the  rate  of  2400  FPM  Is  achieved  from  level  flight  in  only 
5.0  seconds. 

2.4  System  Model 

Figure  2.1  is  a  diagram  of  the  system  model  which  will  be  used  in  the 
following  investigation  of  altitude  tracking.  In  this  diagram  the  Laplacian 
notation  l/s  indicates  the  process  of  time  integration.  The  function  INT  (x) 
is  defined  as  the  value  of  x  truncated  to  the  greatest  integer  less  than  or 
equal  to  x.  Whether  quantization  is  achieved  through  truncation  or  rounding 
is  irrelevant  to  the  rate  tracking  problem.  The  model  employed  here  will 
assume  that  truncation  is  employed.  The  superscript  indicates  a  quantized 
measurement.  The  caret  ("<*")  is  used  to  indicate  an  estimated  quantity.  The 
sampling  rate  is  fixed  at  one  sample  every  t  seconds.  The  nominal  value  of  t 
for  BCAS  is  1  second. 
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3.0  CHARACTERISTICS  OF  QUANTIZED  SAMPLES 

Certain  characteristics  of  quantized  samples  will  now  be  derived  as  a 
prelude  to  the  consideration  of  estimation  techniques.  Mathematical 
relationships  which  are  useful  in  the  following  derivations  are  provided  in 
Appendix  A  of  this  document.  Let  the  altitude  of  an  aircraft  be  given  by 

z(t)  ■  Nq  +  £q  q  +  z  t  (3.1) 


where  N  is  an  integer  and  0  <  eg  <  1.  The  quantity  eg  defines  the  initial 
position  of  the  aircraft  within  the  quantization  level  of  width  q. 

For  simplicity,  the  following  derivation  will  assume  that  z  >  0.  The 
derivations  retain  their  generality  since  the  case  of  z  <  0  is  obtained  by 
merely  reflecting  the  equations  about  the  origin  of  the  coordinate  system.  The 
quantized  altitude  at  the  nc^  sample  instant  can  now  be  written: 

• 

nzt 

z*(nT)  =  Nq  +  q  INT  (eg  +  -  ) 

q 

nx 

«  Nq  +  q  INT  (eg  +  -  )  (3.2) 

T 

where  T  is  the  level  occupancy  time,  i.e.,  the  time  required  for  z(t)  to 
change  by  an  amount  q.  It  can  be  seen  that  the  k11*1  altitude  transition  occurs 
on  sample  n^  which  satisfies  the  following  inequality: 


where  R  is  defined  as  the  function  which  yields  the  fractional  remainder  of 
the  argument,  i.e., 

R(x)  -  x  -  INT(x) 
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The  observed  time  within  a  single  quantization  level,  Tk,  is  related  to 
the  nk  as  follows: 


Tk  “ 


T(nk+i  -  nk) 


k  -  1,2, 3, 4,. . . 


(3.5) 


Note  that  at  least  two  transitions  must  occur  to  produce  the  first 
observation  of  level  occupancy.  Substituting  from  equation  3.4  yields  the 
fact  that 


t  INT  (T/t)  for  R  [(k-eo> - 1  <  1  -  R  (T/t) 


(3.6) 


t  INT  (T/t)  +  t  for  R  [(k-eo> - ]  >  1  -  R  (T/t) 

T 


If  T  is  an  integer  multiple  of  t,  then  the  observed  level  occupancy  time 
will  always  be  equal  to  T.  Otherwise  there  are  two  possible  values  of  Tk*. 
These  two  values  differ  by  t  and  bracket  the  value  of  T.  For  randomly 
selected  initial  conditions  the  shorter  value  occurs  with  probability 
1-R(T/t ).  The  longer  value  occurs  with  probability  R(T/t).  The  value  of  Tk 
does  not  vary  by  more  than  x  for  constant  rate  trajectories.  This  fact  can  be 
used  in  the  detection  of  accelerations  (as  will  be  demonstrated  later). 


4.0  LINEAR  RECURSIVE  TRACKING 

The  original  BCAS  tracking  algorithm  uses  the  alpha-beta  smoothing 
smoothing  equations  which  follow: 


*  z  +  t  z 


n+l,n 


(4.1) 


zn+l,n  +  °  ^zn+l  “  zn+l,n^ 


(4.2) 


zn  +  6 


(«*  ,  -  *  ) 
n+1  n+1 , n 


(4.3) 


r 1 


where  zn+i»n  Indicates  the  value  of  zn+j  as  projected  from  the  time  of 
measurement  zn*.  Once  a  track  Is  firmly  established,  the  values  of  a  and  0 
are  constant  and  these  equations  comprise  a  linear  recursive  tracking 
technique. 


4. 1  Step  Response  of  the  Alpha-beta  Tracker 

The  response  of  the  alpha-beta  tracker  to  a  change  in  altitude  input  of 
one  quantization  level  will  now  be  discussed.  Consider  a  case  in  which  an 
aircraft  with  negligible  altitude  rate  crosses  a  quantization  boundary. 
Assuming  that  the  tracker  estimate  had  previously  converged  to  level  flight  at 
the  reported  altitude,  the  tracker  is  presented  with  a  sudden  discrepancy  of 
magnitude  q  between  the  predicted  and  measured  position.  The  characteristic 
response  is  sketched  in  Figure  4.1.  The  resulting  rate  estimates  can  be 
computed  in  closed  form  as  shown  in  Table  4.1.  For  typical  a  and  $  values, 
the  maximum  rate  error  occurs  on  the  second  or  third  sample  following  the 
transition. 

The  maximum  rate  error  for  1  second  update  Interval  is  plotted  in  Figure 

4.2  for  a  range  of  a  and  0  values.  Normally  choices  of  a  and  0  would  be 
matched  according  to  the  formula 


6  - - 

2-o  (4.4) 

as  suggested  by  Benedict  and  Bordner*.  Points  corresponding  to  this  formula 
are  indicated  in  the  figure. 


4.2  Tracking  Cycle  Behavior 

For  aircraft  climbing  at  low  to  moderate  rates,  several  samples  are 
obtained  at  each  quantization  level.  The  rate  estimates  of  the  recursive 
alpha-beta  tracker  then  tend  to  follow  a  tracking  cycle  in  which  the  velocity 
is  overestimated  at  the  scan  or  scans  immediately  following  a  transition  and 
underestimated  for  later  scans.  Depending  on  the  value  of  the  smoothing 
constants  and  the  rate,  the  estimate  near  the  end  of  the  tracking  cycle  may 
converge  to  zero  or  even  be  opposite  in  sign  to  the  direction  of  the 
transitions.  This  behavior  is  shown  in  Figure  4.3  for  a  climb  rate  of  800 
FPM,  B  -  0. 1. 


♦Benedict,  T.R.  and  Bordner,  G.W. ,  "Synthesis  of  an  Optimal  Set  of  Radar 
Track-While-Scan  Smoothing  Equations,"  IRE  Transactions  on  Automatic 
Control  (July  1962). 
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TABLE  4.1 

RESPONSE  OF  THE  ALPHA-BETA  TRACKER  TO  A  SINGLE  ALTITUDE  TRANSITION 


Scan 

(Altitude  transition 
on  scan  N) 

Position  Response 
After  Update 

Rate  Response 

After  Update 

N 

aq 

Bq 

T 

N  +  1 

(2a  -  a2  -  ag  +  g)q 

Bq 

-  (2  -  a-  B) 

T 

N  +  2 

(3a-  3a2  -5aB  +a3  -  B2 

Bq 

-  (3  _  3a_  4g+  2«e  +  a2  +  B2) 

+  3(5  +  2a2B  +a  B2)q 

T 
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Fig.  4.2.  Maximum  rate  error  of  an  alpha-beta  tracker  due  to  a  single  100  foot 
quantization  level  transition  occuring  in  near-level  flight.  The  errors  shown 
are  proportional  to  q/T  and  hence  the  figure  can  be  adapted  to  any  update  rate 
or  quantization  fineness  by  simply  re-labeling  the  ordinate. 
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Fig.  4.3.  Rate  estimates  of  an  a  -  6  tracker  for  800  FPM  ramp  (a®  0.4  6*0.1) 


5.0  LEVEL  OCCUPANCY  TRACKING 


5.1  Motivation  for  Level  Occupancy  Tracking 

The  underlying  reason  for  the  undesirable  tracker  behaviors  described 
above  is  that  the  error  characteristics  of  the  measurement  system  diverge  from 
those  characteristics  for  which  the  alpha-beta  tracker  is  best  suited.  The 
alpha-beta  tracking  treats  each  data  point  as  an  independent  measurement  (as 
if  errors  were  uncorrelated  from  sample  to  sample).  In  reality,  the  errors 
can  be  highly  correlated,  especially  when  the  rate  is  such  that  multiple 
samples  are  obtained  in  the  process  of  crossing  a  single  quantization  level. 
A  little  thought  reveals  that  all  the  rate  information  is  contained  in  the 
history  of  altitude  transitions.  Increasing  the  sampling  rate  benefits 
tracking  accuracy  not  so  much  by  providing  more  data  points  at  each  level  as 
by  decreasing  the  error  in  the  determination  of  the  level  transition  times. 

The  above  observations  suggest  a  change  in  the  approach  to  altitude 
tracking.  In  order  to  avoid  the  problems  of  the  recursive  update  at  each 
sample  point,  all  data  points  which  fall  within  the  same  altitude  quantization 
bln  are  treated  as  a  single  observation  of  rate.  The  smoothing  equations  are 
then  written  in  terms  of  level  occupancy  time,  T,  the  amount  of  time  required 
for  the  aircraft  to  cross  a  single  quantization  level.  The  true  value  of  this 
time  is 


T 


(5.1) 


The  result  of  the  measurement  process  is  then  viewed  as  a  series  of 
observations  of  past  occupancy  times  plus  an  observation  of  the  time  at  the 
current  level. 

The  difference  between  any  two  successive  altitude  transitions  serves  as 
a  measurement  of  T.  Errors  in  this  measurement  are  attributable  to  the  finite 
sampling  rate.  The  effect  of  these  errors  on  the  rate  estimate  are  reduced  by 
smoothing  successive  values  of  the  level  occupancy  time. 

Note  that  in  this  formulation  of  the  problem,  an  altitude  measurement 
which  yields  the  same  altitude  value  as  the  previous  measurement  does  not 
provide  a  new  measurement  of  T  and  hence  is  not  smoothed.  The  tracking  cycle 
behavior  is  thus  eliminated.  However  the  lack  of  a  transition  may  be 
significant  if  the  aircraft  has  occupied  the  current  level  for  longer  than 
expected  -  in  that  case  the  lack  of  a  transition  may  indicate  that  an 
acceleration  has  occurred  which  requires  that  the  magnitude  of  the  rate 
estimate  be  decreased.  Hence  an  altitude  measurement  without  transition 
requires  a  special  check  to  determine  if  a  transition  is  overdue.  The 
differing  update  procedures  are  Indicated  schematically  in  Figure  5.1. 
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Fig.  5.1.  Basic  update  procedure  for  level  occupancy  tracking 


5.2  Analysis  of  Rate  Tracking  Accuracy 

The  effect  of  quantization  upon  tracking  accuracy  will  now  be  examined  for 
an  algorithm  which  functions  by  smoothing  level  occupancy  time.  The  observed 
value  of  the  level  occupancy  time  is  the  time  difference  between  observed 
altitude  transitions.  Hence  the  observed  occupancy  time  is  always  a  multiple  of 
the  fundamental  sample  rate  t. 

Let  the  first  observed  transition  occur  at  time  t0  when  the  aircraft  Is  a 
distance  eg  q  from  the  level  boundary  most  recently  crossed.  This  situation  Is 
shown  In  Figure  5.2.  Note  that  In  this  notation  0  <  eg  <  t/T  <  1.  Since  the 
direction  of  the  transition  determines  the  sign  of  the  rate,  the  sign  can  be 
viewed  as  a  known  quantity.  Hence  we  may,  without  loss  of  generality,  discuss 
only  the  case  for  z  >  0.  The  position  change  between  the  first  and  second 
transitions  of  the  sampled  data  is 

.  * 

z  Tj  -  q  +  Ejq  -  eQq  (5.2) 


and  hence 

* 

Ti  "  T  (  1  +  el  "  e0>  (5.3) 

Furthermore  It  Is  obvious  that  In  the  general  case 


Tj  “  T  (1  +  Ej  -  ej.j) 


(5.4) 


Consider  the  effect  of  averaging  the  level  occupancy  time  over  k 
observations.  The  estimate  of  T  which  results  can  be  written  as 


I  T 

.  j-LiIj  . 


T  +  T 


(ek  -  eo> 


(5.5) 
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ALTITUDE 


SAMPLED  ALTITUDE  VALUE8 


Notation  employed  to  desc 


Fig.  5.2. 


ribe  a  constant  rate  altitude  profile. 


The  error  in  T  is 


(5.6) 


e 


T 


This  error  is  proportional  to  the  difference  in  the  initial  and  final  values 
of  the  quantization  residual  e  and  is  Inversely  proportional  to  the  number  of 
level  occupancy  transitions  observed. 


The  value  of  e  can  be  no  more  than  t/T  on  the  scan  at  which  a  transition 
is  observed.  Thus  0  <  Ej  <  t/T,  j  ■  0,1,2,...,  and  the  magnitude  of  the  error 
is  limited  according  to 


(5.7) 


Under  the  assumption  that  and  eg  are  independent  and  uniformly 

distributed  over  [0,t/T] ,  the  quantity  e  -  e^  -  Eg  has  the  probability  density 
function 


fe(x) 


T  „  T  T  T 

( - >2  ( - |x|)  - <  X  * - 

T  T  T  T 

(5.8) 

0  Otherwise 


Using  5.8  it  can  readily  be  shown 
is  unbiased.  The  varience  in  T  can  be 


that  E(ex)  ■  0, 
shown  to  be: 


i.e. , 


the  estimate  of  T 


and  hence 


(5.9) 


(5.10) 
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This  error,  expressed  as  a  fraction  of  x  ,  is  Inversely  proportional  to 
the  number  of  observed  level  occupancy  Intervals.  This  expression  is  plotted 
in  Figure  5.3.  It  should  be  noted  that  no  estimation  of  T  is  possible  until 
at  least  two  altitude  transitions  have  been  observed.  The  error  is  never 
greater  than  the  sample  interval  x. 


It  should  be  noted  that  the  error  converges  faster  than  would  be  expected 
for  a  case  where  measurement  errors  were  normally  distributed  white 
noise.  In  the  latter  case,  decreases  according  to  1//  k  rather  than  1/k. 

A  A 

•  ^  • 

If  the  rate  estimate  is  simply  z  *  q/T,  then  the  error  in  2  at  the  kth 
observed  occupancy  time  can  be  expressed  as  a  fraction  of  z  as  follows: 


Z  <tk)  -  2 


k  +  Ek  ~  e  0 


(5.11) 


The  value  of  -  eg  satisfies  the  following  inequality 


x  x 

-  <  ek  "  e0  < 


(5.12) 


Hence  the  maximum  value  of  e«  is: 


k  -  x/T 


(5.13) 


kT  -  1 
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STANDARD  DEVIATION  OF  ESTIMATE, 


NUMBER  OF  OBSERVED  OCCUPANCY  TIMES.  K 


Fig.  5.3.  Summary  of  errors  in  the  estimation  of  T  after  K  observed 
level  occupancy  times.  Note  that  the  number  of  observed  level 
transitions  is  K+l. 
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5.3  Position  Correction 

As  previously  discussed  in  Section  2.2,  it  is  generally  not  possible  to 
determine  exactly  where  an  aircraft  is  within  a  quantization  level.  However, 
it  can  be  seen  from  Figure  5.2  that  on  the  scan  upon  which  a  transition  occurs, 
the  aircraft  must  be  within  a  distance  of  qx  /T  of  the  quantization  boundary 
most  recently  crossed.  Thus,  if  the  boundary  crossed  is  at 
Nq,  (N  -  1,2,3  ...),  then  the  altitude  is  between  Nq  and  Nq  +  z  x .  A  tracking 
algorithm  can  use  this  fact  by  always  correcting  the  position  estimate  so  that 
it  lies  within  this  interval.  In  the  algorithm  described  later,  the  midpoint 
of  the  indicated  altitude  interval  is  selected  as  the  new  estimated  position, 
i.e. , 


z  «  Nq  +  1/2  z  x 

This  correction  is  of  greatest  benefit  when  T  »  x  and  is  of  negligible  benefit 
when  T  =  x  . 

6.0  ALGORITHMIC  PROVISIONS  FOR  TRACKING  RATE  CHANGES 

The  preceeding  section  introduced  the  concept  of  rate  estimation  through 
smoothing  of  level  occupancy  times.  The  results  obtained  by  simple  averaging 
of  observed  occupancy  times  were  analyzed.  Although  this  averaging  technique 
can  form  the  basis  for  estimation,  an  algorithm  is  not  complete  without 
provisions  for  adequate  tracking  in  the  case  of  non-constant  rates.  Three  such 
provisions  are  discussed  in  this  section.  The  first  is  a  provision  for 
decreasing  the  magnitude  of  the  estimated  rate  when  expected  level  transitions 
fail  to  occur.  The  second  is  a  technique  for  preventing  the  gain  of  the 
tracker  from  converging  to  zero  during  long  periods  of  consistent  transitions. 
Finally,  a  procedure  for  testing  the  consistency  of  the  estimated  and  observed 
occupancy  times  is  described.  A  tracking  algorithm  which  combines  all  these 
techniques  is  provi^xd  in  Appendix  B. 

6.1  Update  Procedure  for  Decreases  in  Rate  Magnitude 

As  described  in  Section  5.1,  the  algorithm  under  consideration  smooths  the 
level  occupancy  time  only  when  altitude  transitions  occur.  A  distinct  update 
procedure  must  be  invoked  whenever  expected  transitions  fail  to  occur.  The 
procedure  for  detection  of  overly  long  level  occupancies  is  based  upon  the 
following  test  for  consistency  between  observed  and  estimated  level  occupancy 
times. 

Let  Sj  be  the  difference  between  the  estimated  and  the  observed  value  of  T  on 
the  jth  transition.  That  is 

8j  -  Tj.!  -  T*  j  -  2,3,4....  (6.1) 
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It  has  been  shown  In  equation  3.6  that  the  observed  value  of  T  always 
falls  between  tINT(T/t)  and  tINT(T/t  )  +t.  Since  Tj_i  is  aerely  an 
average  of  previous  observed  level  occupancy  tines  it  is  within  the  sane  range 
of  width  t.  Thus  the  residual  sj  has  magnitude  of  t  or  less: 

I  8 j I  <  T  (6.2) 

Barring  acceleration,  an  observed  occupancy  time  should  never  differ  from 
the  previous  estimate  by  more  than  t.  Hence  if  the  current  time  of  level 
occupancy  (i.e.,  the  time  since  the  last  transition)  is  equal  to  or  greater 

A 

than  T  +  r  ,  then  an  acceleration  which  has  decreased  the  rate  magnitude  is 
indicated.  The  proper  response  of  the  filter  to  this  situation  depends  upon 
the  assumed  aircraft  trajectory  statistics.  A  minimal  response  would  be  to 
decrease  the  rate  estimate  only  as  much  as  required  to  make  it  consistent  with 
the  observed  occupancy  time.  But  it  has  been  found  that  this  results  in  a 
very  slow  convergence  to  zero  which  is  undesirable  when  the  aircraft  has 
leveled  out  at  the  end  of  a  climb  or  descent. 

Under  the  assumption  that  a  return  to  level  flight  is  more  likely  than  a 
transition  to  another  non-level  altitude  rate,  a  better  approach  is  to  force  a 
rather  fast  convergence  to  zero  when  the  data  clearly  indicates  that  an 
acceleration  toward  zero  rate  has  occurred.  This  is  the  approach  taken  in  the 
algorithm  described  in  Appendix  B.  A  variable  which  indicates  excess 
occupancy  time  is  defined  by 


t  (6.3) 

where  T*  is  observed  time  in  the  current  quantization  level  and  T  is  the 
current  estimate  of  the  level  occupancy  time.  Whenever  6  is  greater  than  1, 
the  rate  is  adjusted  toward  zero  according  to  an  experimentally  derived 
formula-  When  &  is  greater  than  5,  it  is  assumed  that  the  aircraft  has 
returned  to  level  flight. 


««  a  -  mt&i  y.iimaaMiHr.* 


6.2  Lower  Limit  on  Tracker  Gain 

The  process  of  averaging  level  occupancy  times  can  be  carried  out  In 
recursive  fashion  by  use  of  the  following  smoothing  equation: 


where 


and 


*»  -  Vi +  s» 


8n  -  1/n 


Tq  -  0 


(6.4) 


(6.5) 


(6.6) 


Note  that  as  the  number  of  smoothing  updates  increases,  the  gain  of  the 
tracker  (as  represented  by  0n)  decreases  toward  zero.  This  is  appropriate 
only  If  the  rate  is  constant.  In  actual  operation,  the  rate  is  never  exactly 
constant  and  hence  the  convergence  of  0  n  to  zero  must  be  slowed  or 
Interrupted. 

A  classic  least  squared  error  approach  to  tracking  in  the  presence  of 
state  perturbations  results  in  a  tracker  gain  which  decreases  toward  a  limit 
which  is  determined  by  the  ratio  of  the  measurement  error  to  the  magnitude  of 
the  perturbation.  In  the  current  case  the  perturbation  Is  related  to  the 
amount  of  velocity  change  which  is  expected  between  smoothing  Instants.  These 
instants  occur  at  time  intervals  of  approximately  T.  There  Is  a  greater 
expected  perturbation  when  T  is  large.  This  implies  that  the  lower  limit  on 
0n  should  be  smaller  for  smaller  T,  larger  for  larger  T.  (Another  way  of 
looking  at  this  is  that  in  order  to  implement  a  "fading  memory"  tracker,  one 
must  assign  less  weight  to  observations  the  further  they  are  in  the  past.  The 
smaller  the  value  of  T,  the  more  recent  were  past  observations,  and  the  more 
heavily  they  should  be  weighted  relative  to  the  current  observation). 
Although  the  true  value  of  T  is  not  available  to  the  tracker,  it  is  sufficient 
to  select  the  limit  on  the  basis  of  T.  In  the  tracking  algorithm  tested  in 
simulation,  a  0  limit  was  selected  according  to  the  following  expression: 
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max 


0.08 


(6.7) 


®min 


_4i_z_nL_ 

T2  +  64 


This  expression  was  developed  by  first  selecting  a  form  which  satisfied 
reasonable  "end  point"  criteria.  Parameter  values  were  then  optimized  through 
simulation.  The  limit  it  imposes  is  plotted  in  Figure  6.1. 


6.3  Consistency  Tests  Applied  to  Level  Transitions 

This  section  discusses  criteria  used  to  determine  whether  the  time  of  an 
observed  level  transition  is  consistent  with  the  existing  rate  estimate. 
Inconsistencies  indicate  that  acceleration  has  occurred  and  that  the  existing 
estimate  may  have  substantial  error.  The  response  to  this  situation  is  either 
to  reinitialize  the  track  or  to  Increase  the  tracker  gain. 


6.3.1  Sign  Consistency 

The  direction  of  the  altitude  transition  is  always  the  same  as  the 
direction  of  the  true  rate.  Note  that  sign  consistency  Is  not  guaranteed  in 
the  simple  alpha-beta  tracker  and  may  result  l)  when  rate  reversals  occur  or 
2)  as  a  consequence  of  the  tracking  cycle  behavior  discussed  in  4.2  as  the 
rate  estimate  oscillates  around  zero.  In  the  algorithm  formulated  in  Appendix 
B,  sign  consistency  is  required  at  all  times.  This  means  that  whenever  the 
sign  of  the  transition  is  opposite  to  the  sign  of  the  rate  estimate,  the  track 
is  reinitialized. 
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Fig.  6.1.  Lower  limit  on  smoothing  gain. 
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6.3.2  Single-scan  Consistency  Test 


As  shown  in  Section  6.1,  the  magnitude  of  the  residual 


(6.8) 


should  never  exceed  t.  Hence  If  |sjl  is  greater  than  t,  acceleration  has 
occurred  and  reinitialization  is  desirable.  When  such  reinitialization 
occurs,  averaging  begins  anew  and  smoothing  parameter  Dr  of  Equation  6.4 
should  be  reset  to  the  value  for  n*l. 

6.3.3  Test  of  Summed  Residuals 

If  an  acceleration  alters  the  value  of  T  by  an  amount  less  than  T,  the 
residual  may  never  exceed  x  (or  may  not  exceed  t  until  several  level  occupancy 
intervals  have  passed).  The  single  scan  consistency  test  will  not  result  tn 
detection.  But  acceleration  may  still  be  detectable  by  a  test  which  sums 
residuals  over  more  than  one  transition.  This  test  is  based^upon  the  fact  that 
the  sum  of  the  residuals  tends  toward  zero  except  when  T  is  in  error.  The 
test  is  implemented  by  computing,  at  each  smoothing,  a  weighted  sum  of  the 
residuals  according  to  the  formula 


8j  -  If  8j_j  +  sj  (6.9) 

Whenever  Ij  exceeds  a  certain  threshold,  an  excess  residual  is  declared. 
The  value  of  i<  is  reset  once  such  a  detection  has  occurred.  The  parameter  y 
is  set  to  slightly  less  than  unity  to  provide  a  gradual  reset  of  Sj  in  periods 
within  which  no  detection  occurs.  In  the  algorithm  described  in 
Appendix  B,  the  value  of  y  is  0.8  and  the  detection  threshold  for  Ij  is  1.35. 
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7.0  SIMULATION  RESULTS  FOR  A  LEVEL  OCCUPANCY  TRACKER 


The  performance  benefits  to  be  derived  from  the  tracking  techniques 
developed  In  the  previous  sections  cannot  be  adequately  evaluated  without 
testing  a  complete  algorithm  which  properly  integrates  the  various  tracking 
features.  This  section  Introduces  one  such  integrated  algorithm  and  presents 
simulation  results  which  demonstrate  its  performance  in  a  variety  of 
situations. 


7.1  Description  of  the  Tracking  Algorithm 

A  complete  description  of  the  algorithm  used  for  simulation  is  provided 
in  Appendix  B.  Several  details  of  this  algorithm  (including  the  choice  of 
parameter  values)  were  derived  by  experimentation  with  a  range  of  options.  No 
attempt  is  made  here  to  present  the  intermediate  results  which  led  to  this 
final  form.  However,  careful  inspection  of  the  simulation  data  supports  the 
contention  that  the  final  algorithmic  form  is  "near  optimum"  in  that  It 
achieves  almost  all  the  performance  improvement  which  can  be  expected  from 
level  occupancy  tracking. 

The  basic  features  of  the  algorithm  can  be  summarized  as  follows: 

a.  A  single  level  transition  following  an  extended  period  of 
constant  altitude  flight  results  in  initialization  of  the  rate 
estimate  magnitude  to  a  nominal  value  of  480  FPM  (note  parameter 
Pi).  This  rate  decays  by  10%  (note  parameter  P3)  on  each 
successive  scan  without  a  transition. 

b.  The  second  transition,  if  consistent  in  sign  with  the  first, 
results  in  initialization  of  T  to  the  observed  time  between 
transitions.  No  routine  decay  is  then  allowed. 

c.  All  subsequent  transitions  which  satisfy  consistency  tests 
result  in  recursive  averaging  of  the  observed  bin  occupancy 
times. 

d.  If  the  time  for  which  a  level  has  been  occupied  exceeds  by  1.5x 
to  5t  (note  parameters  P5  and  P6)  the  time  predicted  by  T» 
the  rate  is  driven  toward  zero  by  an  empirically  determined 
formula.  If  the  excess  time  is  greater  than  5t  ,  a  return  to 
constant  altitude  flight  (zero  rate)  is  effected. 

e.  If  at  any  point  an  inconsistency  is  noted  in  the  sign  of  the 
rate  estimate  and  the  direction  of  a  transition,  the  rate 
estimate  is  re-initialized  as  in  (a). 
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f.  If  at  the  time  of  a  transition,  the  observed  occupancy  time 

differs  from  the  estimate  by  more  than  t,  the  rate  is 

reinitialized  to  the  rate  corresponding  to  the  observed 

occupancy  time.  This  provides  a  quick  response  to  a  major 
change  in  an  established  rate. 

g.  If  the  smoothed  residual  of  equation  6.9  exceeds  t  at  the  time 
of  any  transition,  the  gain  of  the  tracker  is  increased  and  the 
convergence  of  the  gain  towards  its  lower  limit  is  set  back. 

7.2  Simulation  Results  at  l-Second  Update  Interval 


The  series  of  figures  which  follow  compare  the  rate  tracking  performance 
of  the  level  occupancy  tracker  with  simple  alpha-beta  tracking.  The 
simulation  employed  no  measurement  jitter  and  used  a  constant  update  Interval 
of  1  second.  All  aircraft  rate  changes  took  place  with  an  acceleration 
magnitude  of  0.25g  (8ft/sec^). 


7.2.1  Step  Response 

The  response  of  the  trackers  to  a  single  Isolated  altitude  transition  (a 
"step"  in  altitude)  is  shown  in  Figure  7.1.  As  previously  noted,  the  step 
response  of  the  level  occupancy  tracker  is  arbitrary,  being  determined  by  the 
parameters  Pi  and  P3.  In  BCAS  testing  the  response  of  the  B  -  0. 1  tracker  was 
found  to  be  unacceptable  In  this  situation.  A  value  of  B  *  0.05  was  later 
used  for  BCAS  alpha-beta  tracking  and  found  to  offer  acceptable  step  function 
response.  It  can  be  seen  that  the  level  occupancy  tracker  appears  to  be  of 
essentially  equivalent  acceptability. 


7.2.2  Ramp  Response 

The  response  of  the  trackers  to  various  simple  climb  profiles  will  now  be 
shown.  A  "ramp"  climb  profile  will  be  defined  as  a  profile  which  involves  an 
initial  period  of  level  flight,  acceleration  to  a  specified  rate,  a  period  of 
constant  rate  climb,  acceleration  back  to  zero  rate,  and  a  final  period  of 
level  flight.  Generally  the  largest  tracking  errors  will  occur  during  or 
immediately  following  the  periods  of  acceleration.  The  magnitude  of  the 
errors  is  strongly  dependent  upon  the  ramp  rate  magnitude.  Figure  7.2  plots 
rate  estimates  for  a  450  FPM  ramp.  The  tracking  cycle  behavior  of  the 
alpha-beta  tracker  is  quite  obvious  here,  even  with  the  reduced  tracker  gain. 
Note  that  the  level  occupancy  tracker  eliminates  the  tracking  cycle  behavior 
after  the  second  transition.  When  the  climb  terminates,  the  drifting  of  the 
alpha-beta  tracker  back  toward  zero  results  in  a  fortuitous  convergence  toward 
the  true  rate. 


Figure  7.3  plots  the  ramp  response  for  an  800  FPM  ramp.  Here  the 
alpha-beta  parameters  are  well  matched  to  the  rate  and  convergence  for  both 
trackers  seems  to  occur  at  about  the  second  transition.  The  tracking  cycle 
behavior  is  noticeable,  but  is  a  smaller  fraction  of  the  true  rate. 
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Fig.  7.2.  Rate  estimation  performance  for  a  450  FPM  ramp  profile. 
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Fig.  7.3.  Rate  estimation  performance  for  a  800  FPM  ramp  profile. 
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Figure  7.4  provides  results  for  a  2100  FPM  ramp.  At  this  higher  rate  it 
can  be  seen  that  the  alpha-beta  tracker  does  not  converge  to  the  true  rate 
until  about  the  fifth  altitude  transition  (13  seconds  after  the  first 
transition).  There  is  a  period  of  time  after  the  second  transition  when  the 
rate  is  severely  underestimated  in  comparison  to  the  level  occupancy  tracker. 
Similarly,  the  alpha-beta  tracker  requires  many  scans  to  recover  when  the 
climb  terminates. 

Figure  7.5  presents  results  for  a  5000  FPM  ramp.  For  this  high  final 
rate  the  period  of  acceleration  is  prolonged  and  several  transitions  occur 
during  the  acceleration.  The  level  occupancy  tracker  appears  to  do  better 
than  the  alpha-beta  tracker  during  and  following  the  acceleration  periods. 

One  meaningful  measure  of  performance  is  the  number  of  scans  on  which  the 
magnitude  of  the  rate  error  exceeds  a  given  threshold  of  significance.  This 
measure  indicates  the  amount  of  time  for  which  the  system  is  vulnerable  to 
failures  caused  by  tracking  error.  Figure  7.6  summarizes  the  ramp  performance 
of  the  alpha-beta  and  level  occupancy  trackers  using  such  an  error  count.  The 
number  of  scans  on  which  the  magnitude  of  the  rate  estimation  error  exceeded 
600  FPM  is  given  for  a  range  of  ramp  rates.  For  rates  below  1000  FPM,  this 
measure  shows  little  apparent  differences  in  tracker  performance.  But  as 
rates  exceed  1500  FPM,  the  advantage  of  the  level  occupancy  tracker  becomes 
significant. 

It  should  be  noted  that  in  most  cases  a  change  in  reported  altitude  of 
only  200  feet  (2  transitions)  is  required  in  order  for  the  tracker  to  converge 
to  a  reasonable  estimate  of  the  vertical  rate.  Hence  for  the  level  occupancy 
tracker  it  is  more  appropriate  to  think  in  terms  of  the  number  of  transitions 
required  for  convergence  rather  than  the  length  of  time  required.  The  tracker 
converges  very  rapidly  when  the  vertical  rate  is  high,  since  high  rates 
quickly  produce  the  required  number  of  altitude  level  transitions. 

7.2.3  Steady  State  Performance 

RMS  rate  errors  for  steady  state  climbs  at  various  rates  are  shown  in 
Figure  7.7.  Although  the  steady  state  errors  are  never  large  for  either 
tracker,  the  level  occupancy  tracker  exhibits  smaller  errors,  especially  for 
lower  rates  (below  1500  FPM). 

8.0  TRACKING  WITH  LESS  FREQUENTLY  SAMPLED  DATA 

Ground-based  air  traffic  control  sensors  typically  possess  update 
intervals  of  more  than  4  seconds.  At  these  sampling  rates,  aircraft  with  rate 
magnitudes  of  more  than  1500  feet  per  minute  will  cross  more  than  one 
quantization  level  between  samples.  Several  modifications  must  be  made  to  the 
tracking  algorithm  described  in  Appendix  B  in  order  to  apply  it  to  rate 
tracking  in  this  regime: 

-  The  value  of  the  algorithmic  parameters  must  be  adjusted. 

-  The  smoothing  equation  must  be  generalized  to  accomodate  transitions  of 
more  than  one  quantization  level. 

-  The  consistency  testa  must  be  extended  to  consider  the  number  of 
transitions  between  samples  as  well  as  the  number  of  samples  between 
transitions. 
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Fig.  7.4.  Rate  estimation  performance  for  a  2100  FPM  ramp  profile. 
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Fig.  7.5.  Rate  estimation  performance  for  a  5000  FPM  ramp  profile. 
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An  algorithm  incorporating  these  changes  is  provided  in  Appendix  C.  The 
principal  logic  changes  incorporated  into  this  algorithm  will  be  described  in 
the  next  two  subsections.  Simulation  results  for  a  4.7  second  data  rate 
(typical  for  terminal  air  traffic  control  radars)  will  then  be  presented. 

8.1  Algorithm  Modifications 

8.1.1  Generalized  Smoothing  Equation 

The  update  procedure  for  estimation  of  the  level  occupancy  time  can  be 
generalized  to  handle  cases  in  which  more  than  one  quantization  level  is 
crossed  between  samples.  The  basic  approach  is  still  to  average  the  level 
occupancy  times.  If  M  transitions  have  occurred  in  a  single  sample  interval 
of  duration  t,  then  the  tracker  will  respond  as  if  M  single  transitions  were 

Recall  that  in  the  averaging  process 

n-I 

I  T*  (8.1) 

j-l__ 

n-1 

and  for  the  current  update,  equation  (6.4)  is  employed  with  6  “  1/n.  Hence 
for  a  transition  across  M  levels,  the  average  after  update  should  reflect 
n-l4M  transitions  in  a  time  Interval  which  has  been  incremented  by  t  over  the 
time  of  the  previous  update.  Using  equation  (6.4)  this  yields 


■  -rrrmnr  +  6  <T 

Note  that  in  the  case  of  M*>1,  this  formula  is  equivalent  to  equation  6.4. 

8.1.2  Extended  Consistency  Test 
Using  the  altitude  expression  provided  in  equation  3.1,  i.e. 
z(t)  •  N  q  +  Eg  q  +  z  t, 

the  number  of  transitions  occuring  between  sample  n-1  and  sample  n  can  be 
written 


M>  1 


M>  1  (8.2) 


reported  at  intervals  of  t/M. 


n-1 


n  x  n  x  -  x 

Mn  -  INT  e0  + -  -  INT  eQ  + - 

T  T 
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Using  Theorem  5  of  Appendix  A,  this  expression  can  be  written 


^n  " 


T 


INT  ( - ) 

T 


when 


,  n  T. 

R(e  0  + - )  <  1 

T 


(8.3) 


INT  (---) 


T 


+  1 


when  R(e  g  + - )  > 

T 


A  consistency  test  can  be  based  upon  in  a  manner  analagous  to  the 
consistency  test  based  upon  level  occupancy  times  as  described  in  Section  6. 
Define  a  residual  Sj  according  to 


-  Mj 


(8.4) 


This  quantity  is  the  expected  number  of  level  transitions  between  samples  less 
the  number  of  transitions  which  actually  occurred.  In  the  absence  of 
acceleration  this  residual  should  average  to  approximately  zero  and  should 
have  maximum  magnitude  of  1.  The  residual  becomes  more  positive  when  there  is 
an  acceleration  which  decreases  the  rate  magnitude  and  becomes  more  negative 
in  the  presence  of  acceleration  which  increases  the  rate  magnitude.  A  summed 
value  of  the  residual  defined  by 

-  y  Vi  +  sj 

can  be  tested  in  order  to  detect  accelerations.  When  |sjl  exceeds  a  selected 
threshold,  the  tracker  gain  is  increased  to  better  follow  the  apparent 
acceleration. 

8.2  Simulation  Results  at  4.7-Second  Update  Interval 

Appendix  C  contains  a  listing  for  a  tracking  algorithm  modified  to 
function  at  a  4.7-second  update  rate  according  to  the  principles  introduced 
in  Section  8.1.  This  is  a  rate  that  is  typical  of  air  traffic  control  radars. 
Simulation  results  for  several  altitude  profiles  are  provided  in  Figures  8.1 
through  8.4,  The  parameter  values  used  for  the  alpha-beta  tracker  (a  •  0.464, 
8  ■  0.144)  are  those  employed  in  the  proposed  tracker  for  the  Automatic 
Traffic  Advisory  and  Resolution  System  (ATARS).  It  can  be  seen  that  although 
the  step  function  response  and  tracking  cycle  behavior  are  not  considerations 
at  the  longer  update  interval,  the  response  to  acceleration  can  still  be 
Improved  through  use  of  a  level  occupancy  tracker.  Figure  8. 5  provides  a 
comparison  of  the  error  counts  of  the  two  trackers  for  ramp  profiles  at 
various  rates. 
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Fig,  8.1.  Rate  estimation  performance  for  single  altitude  transition, 
4,7  -  second  update  interval. 
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Fig.  8.4.  Rate  estimation  performance  for  5000  FPM  ramp  profile, 
4.7  -  second  update  interval. 
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Fig.  8.5.  Number  of  scans  on  which  the  rate  estimation  error  exceeded 
600  FPM  for  various  ramp  profiles  for  4.7  -  second  update  interval. 


9.0  SUMMARY  AND  CONCLUSIONS 

The  simple  linear  recursive  algorithms  typically  employed  in  the  tracking 
of  aircraft  vertical  motion  exhibit  undesirable  rate  responses  following 
widely  spaced  transitions  between  adjacent  altitude  quantization  levels.  The 
selection  of  tracking  gain  also  involves  a  difficult  trade-off  between  steady 
state  tracking  performance  and  response  to  acceleration.  An  alternative 
approach  to  tracking  is  examined  in  this  document.  Instead  of  smoothing  each 
altitude  report  in  an  identical  manner,  the  alternative  tracker  smooths  the 
time  the  aircraft  spends  at  each  quantization  level.  This  level  occupancy 
tracking  results  in  a  rate  estimate  which  is  not  influenced  by  redundant 
samples  taken  while  the  aircraft  is  within  the  same  quantization  level.  Since 
two  level  transitions  are  necessary  to  measure  level  occupancy  time,  the 
response  to  a  single  Isolated  transition  is  arbitrary  and  is  controlled  to 
eliminate  problems  of  over-response  in  near-level  flight.  Further  Improvement 
in  tracking  performance  results  from  the  use  of  consistency  tests  which  use 
properties  of  quantized  data  to  detect  inconsistencies  between  tracker 
estimates  and  observed  data.  When  inconsistencies  are  discovered,  an 
Immediate  correction  is  applied  to  the  estimate  and  the  tracking  gain  is 
adjusted  accordingly. 

Simulation  of  a  particular  realization  of  this  type  of  tracking  has  been 
conducted  at  l-second  and  4.7-second  update  rates.  The  simulation  indicates 
that  the  alternative  tracking  algorithm  can  be  Independently  optimized  to 
perform  well  in  cases  of  isolated  transitions,  steady  state  and  accelerating 
trajectories.  The  most  significant  performance  difference  between  the 
conventional  alpha-beta  tracking  algorithms  and  the  alternative  algorithm  is 
in  the  response  to  the  Initiation  or  termination  of  altitude  rates  above  1500 
FPM  magnitude. 
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APPENDIX  A 


MATHEMATICAL  RELATIONSHIPS 


This  appendix  presents  certain  mathematical  definitions  and  relationships 
which  are  useful  in  the  analysis  of  quantization  effects.  In  the  following 
equations  the  quantities  m  and  n  are  integers. 

Definition:  INT(x)  *  the  largest  integer  less  than  or  equal  to  x. 

Examples :  INTC3.6)  -  3.0,  INT(3.0)  -  3.0,  INT  (-.6)  -  -  1.0. 


Definition:  R(x)  •  x  -  INT(x) 

Examples :  R(3.6)  *  0.6,  R(3.0)  “0.,  R(-.6)  “  0.4 

~~  Note  that  for  x  >  0,  R(x)  is  the  fractional  part  of  x. 


Theorem  1 : 

0  <  R(x)  <  1 

Theorem  2: 

lNT(m  +  x)  -  m  +  INT(x) 

Theorem  3: 

R(n  +  x)  -  R(x) 

Theorem  4: 

R(n  -  x)  ■  R[l-R(x)] 

Theorem  5: 

INT(x)  when  1  -  R(x)  >  R(y) 

INT(x  +  y)  -  INT(y)  - 

INT(x)  +  1  when  1  -  R(x)  <  R(y) 


Theorem  6: 


R(x  +  y) 
or  equivalently 
R(x  +  y) 


R(x)  +  R(y)  -  INT  [R(x)  +  R(y)] 


R(x)  +  R(y)  when  R(x)  +  R(y)  <  1 

R(x)  +  R(y)  -  1  when  R(x)  +  R(y)  >  1 
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APPENDIX  B 


ALGORITHM  FOR  1-SECOND  UPDATE  INTERVAL 


This  appendix  provides  details  on  the  level  occupancy  tracking  algorithm 
used  to  generate  the  simulation  results  for  a  1.0  second  update  rate  (see 
Section  7.0).  An  overview  of  the  principal  sections  of  the  logic  is  provided 
in  Figure  B.l.  Tables  B. 1  and  B.2  provide  definitions  of  variables  in  the 
track  file  and  parameters  used  in  the  logic.  The  logic  was  divided  into  two 
FORTRAN  subroutines.  The  first  subroutine  (see  Figure  B.2)  initializes  new 
tracks.  In  BCAS  a  track  initiated  by  the  BCAS  surveillance  function  exists 
prior  to  activation  of  the  collision  avoidance  logic  and  hence  this  track  is 
used  for  initialization.  But  the  simple  surveillance  tracker  does  not  store 
all  the  quantities  needed  to  fully  initialize  the  level  occupancy  tracker.  In 
particular,  it  does  not  store  the  time  of  the  preceding  altitude  transition. 
Hence  those  portions  of  the  logic  which  test  level  occupancy  must  be  by-passed 
until  either  the  second  altitude  transition  is  observed  or  enough  time  (about 
18  seconds)  has  passed  to  guarantee  that  the  aircraft  is  in  near  level  flight. 
During  the  interim  period,  the  "start-up"  logic  employs  simple  alpha-beta 
smoothing  equations. 

A  FORTRAN  listing  for  the  update  logic  is  provided  in  Figure  B. 3.  The 
following  notes  will  be  helpful  in  implementation  of  this  algorithm: 

1)  The  quantities  Q  and  DT  are  treated  as  parameters  in  the  software,  but 
they  can  be  replaced  by  their  fixed  values  (Q  ■  100  feet  and  DT  <*  1.0  second). 

2)  As  currently  written,  the  state  vector  elements  ZM0D(8)  and  ZM0D(9)  are 
not  used  simultaneously.  Hence  it  would  be  possible  to  revise  the  code  to 
treat  them  as  the  same  quantity,  thus  reducing  the  size  of  the  aircraft  state 
vector  by  one  element. 

3)  The  quantity  ZMOD(6)  is  the  time  of  last  track  update.  If  the  track  is 
always  updated  at  a  specified  rate,  this  quantity  need  not  be  stored.  This 
would  reduce  the  size  of  the  aircraft  state  vector  by  one  element. 

4)  The  following  quantities  are  internal  to  the  subroutines  and  are 
defined  foe  computational  convenience. 


BETA1 

ISGN 

BLIM 

QSIGN 

DBINS 

TCUR 

DELT 

TEST 

OZM 

TNDEX 

3)  In  the  code  presented  here  an  altitude  report  of  zero  is  used  to 
indicate  missing  data.  Actually,  zero  is  a  valid  Mode-C  report  value.  In 
implementing  this  code,  another  default  value  should  be  chosen. 
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Fig.  B.l.  Overview  of  tracking  algorithm  showing  principal  modules. 
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TABLE  B. 1 


ZMOD(L) 
ZMOD(2) 
ZMOD( 3) 
ZMOD(4) 
ZMOD( 5) 
ZMOD(6) 
ZMOD(7) 
ZMOD(8) 

ZMOD(9) 

ZMOD(IO) 

ZM 


CONTENTS  OF  THE  TRACK  FILE 


Estimated  Altitude  (ft) 

Estimated  altitude  rate  (fps) 

Time  last  Mode-C  report  was  received  (sec) 

Previously  reported  altitude  (ft) 

Time  of  transition  to  previously  reported  altitude  (sec) 

Time  of  last  track  update  (sec) 

Estimated  level  occupancy  time  (sec) 

Firmness  of  rate.  If  equal  to  zero,  Indicates  rate  Is  based 
upon  assumption  of  level  flight  or  observation  of  a  single 
altitude  transition.  If  equal  to  1  or  more,  it 
equals  the  number  of  observed  occupancy  times  for  the 
current  rate  (but  it  may  be  reset  by  consistency  tests). 

Start-up  Counter.  Used  in  establishing  track. 

Summed  residual.  Used  to  detect  a  trend  In  the 

tracker  residuals  (which  indicates  vertical  acceleration). 

Mode-C  altitude  report.  Set  to  0  when  no  report  has  been 
received. 


TABLE  B.2 


PARAMETERS  USED  IN  ALTITUDE  TRACKING 


VARIABLE 

NAME 

DEFINITION 

NOMINAL 

VALUE 

DT 

Nominal  Time  Between  Updates 

1.0  sec.  (BCAS) 

Q 

Quantization  Bin  Width 

100  ft. 

PI 

Magnitude  of  Rate  Allowed  Following 

Isolated  Altitude  Transition 

8  FPS 

P3 

Decay  Factor  When  No  Reinforcing 

Transition  has  Occurred 

0.90 

P4 

Stiff  Rate  Smoothing  Parameter 

0.04 

P5 

Excess  Bin  Occupancy  Time  Which  Results 
in  Transition  to  Level  Flight 

5.0  sec. 

P6 

Excess  Bin  Occupancy  Time  Which  Results 
in  Correction  to  Altitude  Rate  (units  of  DT) 

1.5 

P7 

Amount  of  Discrepancy  in  Bin  Occupancy 

Times  Which  Triggers  Reinitialization  of 

Tracker  Vertical  Rate  (units  of  DT) 

1.5 

P8 

Parameter  Used  to  Position  an  Estimated 

Bin  Transition  Time  Within  an  Interval 
of  Missing  Data 

0.6 

P9 

Position  Smoothing  Parameter 

0.3 

P10 

Smoothing  Gair  Used  To  Compute  Summed 

Residual,  ZMOJ(IO) 

0.80 

P1I 

Value  of  Bin  Occupancy  Smoothing  Parameter 

Used  When  Excess  Residuals  are  Detected 

0.70 

P12 

Reset  Magnitude  for  the  Summed  Residual, 
ZMOD(IO) 

0.30 

P13 

Value  of  ZM0D(9)  at  Which  Transition  From 
Start-up  Smoothing  to  Normal  Smoothing  Occurs 

18 

PI4 

Threshold  Magnitude  for  ZMOD(IO).  Used  to 
Detect  Excess  Summed  Residual 

1.35 
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Fig.  B.2.  Logic  used  to  initialize  a  track  based  upon  a 
previous  existing  track  with  altitude  ZAB(l)  and  altitude 
rate  ZAB(2).  T  is  current  time  and  ZM  is  the  currently 
reported  altitude. 
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SUBROUTINE  VTRMOD  1 

C  -  -  -  -  MODIFIED  VERTICAL  TRACKING  11-19-80  JWA  8 

COMMON /COM A/ ALFA .BETA . T.0T.ZM.ZAC<6).ZAB(6>,  3 

K  RMSI4  )  ,  IOUT  .KOEV  .  10TRAJ  .KT  I  TLE  (  1  0  )  .Dl  AG  (20  I  * 

COMMON/COMMOO/PARM .OTM.Q.ZDINIT. Z MODI  I  8 ). ALFA  I  , BET A  I  ,  5 

K  IPATHC80)  ,KPATH(80)  ,P1  , PS , P3 . P4 , P5 . P6 . P7 . PS .  6 

K  P9.PIO.P1 1 .PI8.P13.P14  7 

C-----------------------  8 

IF  ( T . GT . I . )  GO  TO  808  9 

C--  .  ----------  PARAMETER  SETTING  -------  10 

P I -B .  II 

P 3-0.90  18 

P4- . 04  13 

P5*5.  1 4 

P6-1 .5  15 

P7- 1 .5  16 

P8-.6  17 

P9-.3  18 

PI  0* . 80  19 

P 1 1 > . 7  80 

PIS- . 3  81 

PI3-18.  88 

P I  4- 1  . 35  83 

808  CONTINUE  8H 

C----------------------  UPDATE  -------  25 

ZP-ZM00(l)t(T-ZM0D(6l)«ZM0Dt8)  86 

IF  (ZM.GT.O.)  GO  TO  490  87 

C-  ------  ------  ALTITUDE  REPORT  MISSING:  COAST  TRACK  -  -  -  -  28 

ZM0D11I-ZP  89 

GO  TO  808  30 

V  90  CONTINUE  31 

DZM-ZM-ZMOO I 4 I  38 

DBINS-ABS1DZM1 /Q  33 

ZMOD (9) -ZMOD 19  I ♦ 1  . ♦ 1 0«DB I  NS  34 

IF  ( ZMOD ( 9)  . GT . P I  3)  GO  TO  495  35 

C--. --------------  START-UP  SMOOTHING  ---------  36 

ZMOOI I > »ZP*P9* ( ZM-ZP)  37 

ZMOO(8I*ZMOO(8!*P4«(ZM-ZP)/(T-ZMOD(6))  38 

IF  (OZM.NE.O.l  GO  TO  590  39 

GO  TO  808  40 

495  CONTINUE  41 

I F  (OZM . EQ . 0 .  1  GO  TO  70  I  42 

C-------------------  -  -TRANSITION  LOGIC  ------  43 

I  SON* I  NT ( S I  ON (  I  .  , DZM I)  44 

QSIGN-QMSGN  45 

TEST-ZMOO < 8 ) «DZM  46 

IF  < TEST . LE . 1 00 . I  GO  TO  580  47 

c------------------  TREND  EXISTS  -  -  -  -.-  -  -  48 

TPREV- ( T-ZMODI 5) ) /OBINS  49 

DEL T-TPREV-ZMOD < 7 )  50 


Fig.  B.3.  FORTRAN  subroutine  used  to  update  a  track  at  1  second  intervals. 
This  routine  is  called  once  each  scan  following  the  scan  of  initialization. 
As  a  programing  convenience,  the  setting  of  parameters  has  been  Incorporated 
into  this  sub  routine. 
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IF  I ZMOO (8 ) . LE . 0 . I  GO  TO  540  51 

IF  I ABS (OEL T I . GT . P7»OT )  GO  TO  540  52 

ZMOD  I  IO)«PIO*ZMOD<  1 0  >  *DEL  T  53 

01  AG  I  1) -ZMODl  10)  54 

IF  (ABSIZMOOt 10) ) .LE .P14)  GO  TO  530  55 

C  .  -------------  EXCESS  RESIDUAL  DETECTED-  -  -  56 

BE  T  A l »P 1  I  57 

ZMOO (81*3.  58 

ZMOO< 10)-SIGN(PI2.ZHO0( 10) )  59 

GO  TO  533  60 

530  CONTINUE  61 

C  .  ---------  CURRENT  RATE  WITHIN  LIMITS  :  SMOOTH  -  -  62 

Z7-ZM00I7I  63 

BL I M* ( Z7- 1 . ) • *2/ I Z7* *2*64 . )  64 

BETAI-AMAXIII./I ZMOO  <  8 1 ♦ 1  .  )  . BL I M .  . OB )  65 

ZM00(8> *ZM00(8) ♦ I .  66 

533  CONTINUE  67 

ZM0DC7)-ZM00(7)«-8ETA1*(TPREV-ZM0D(7))  68 

ZMOO C 2 ) -QS I GN/ ZMOO t 7  )  69 

ZMOO I  1  I -ZP«P9* I ZM- ZP )  70 

GO  TO  590  71 

540  CONTINUE  72 

C  -  -  -------  .  -  REINITIALIZE  RATE  ------  73 

ZMOD I  7 )  -AMAX I  (  TPRE V  ,  I  . 4 )  74 

ZMOD l 2 ) -QS I GN/ ZMOO ( 7 )  75 

ZMODllOI-O.  76 

ZMOD (8 ) -  1  .  77 

Z  MOO (  l  ) -ZM-QSIGN/2*ZM00(2) •  DT  /  g  78 

GO  TO  590  79 

580  CONTINUE  BO 

SET  RATE  TO  VALUE  FOR  SINGLE  TRANSITION  -  81 

ZMOO < 2 ) -P 1 • I SGN  82 

ZMOO ( I )-ZM-QSIGN/c*ZM00(5) »DT/2  83 

ZMOD ( 7 ) *QS I GN/ ZMOO ( 2 1  84 

ZMOO ( 8 ) - 0  .  85 

ZMOO!  I  0 ) -0  .  86 

590  CONTINUE  87 

c-  ------------  UPDATE  LEVEL  TRANSITION  OAT*  ------  88 

ZMOO I  4 ) ■ ZM  89 

ZH00(5)-T  90 

IF  ( ZMOD l 3 ) . GE . ZMOO ( 6 ) I  GO  TO  802  91 

C  -  -  -  CORRECT  TRANSITION  TO  FALL  WITHIN  PERIOD  OF  MISSING  OATA  -  -  -  92 

ZMOO I  5 ) *T*P8* I ZMOO I  3  I -T*OT >  93 

GO  TO  802  94 

C . -------- . .  95 

C-. NO  TRANSITION  LOGIC  -------  96 

701  CONTINUE  97 

ZMOOI I ) • ZP+P9* (ZM-ZP)  98 

TCUR»T-ZMOD (5 ) +OT  99 

TNOEX* ( TCUR-ZMOOI 7 ) ) /OT  100 

Fig.  B.3.  Continued. 
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IF  (TN0EX.GT.P5)  GO  TO  630  101 

IF  (TNOEX.GE.P6)  GO  TO  610  102 

C---  --------  NORMAL  UPDATE  -  NO  TRANSITION  EXPECTED  -----  103 

c-  ------  ---------  RATE  DECAYS  ACCORDING  TO  P3  .  104 

IF  ( ZMOD ( 8 1 . GE . 1 . )  GO  TO  802  105 

ZMOD < 2 ) -ZMOO ( 2 > *P3  106 

ZMOD < 7  I «Q/ ( ABS (ZM00(2 I  )  +  .  I  I  107 

GO  TO  002  100 

610  CONTINUE  109 

C  -  -  -  -  SIN  OCCUPANCY  LONGER  THAN  EXPECTEO  -  ALTER  EXTERNAL  RATE  -  -  -  110 

T7-ZM0D(7>*< . 3* ZMOD (7)+.5*0T)*(TNDEX-.3)«»2  111 

ZM0D(2i-SIGN(Q.ZM00(21 1 /T7  112 

ZMOD < 8 ) -AMAX 1  (2  .  . ZM0018 1  -  1  .  1  113 

GO  TO  802  1  I  4 

630  CONT INUE  1  15 

C----. -------------  TRANSITION  TO  LEVEL  FLIGHT  -  -  -  116 

ZMOD ( 1 ) -ZM  117 

ZMOD (21-0.  1  1  0 

ZMOD ( 7 1 -99.  I  19 

ZMOO (81*0.  120 

ZMOO (  1 0  )  -0  .  12  1 

802  CONTINUE  122 

ZMOO < 8) -AMINI ( ZM00(8 I , 1 0 . I  123 

IF  (ZM.GT.O.)  ZMOO ( 3 ) ■ T  124 

ZMOD ( 6 1 • T  125 

RETURN  126 

END  127 

Fig.  B.3.  Continued. 
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APPENDIX  C 


ALGORITHM  FOR  4.7-SECOND  UPDATE  INTERVAL 


A  FORTRAN  subroutine  used  to  implement  a  level  occupancy  tracker  for  a 
4.7  second  update  interval  is  given  in  Figure  C.l.  Variable  definitions  and 
further  explanation  of  the  algorithmic  structure  can  be  found  in  Appendix  B 
and  Section  8.0  of  this  document. 
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JWA  1 

SUBROUTINE  VTRNOD  I 

C  -  -  -  -  MODIFIED  VERTICAL  TRACKING  FOR  4.7  SECONOS  UPOATE  RATE  2 

C0MM0N/C0MA/ALFA.B£TA.T.0T.ZM.ZAC<6I . TAB ( 6 ) ,  3 

K  RMSC4 ) . IOUT ,KOEV , IOTRAJ. KT I TLE (I 0 ) . Dl AO ( 20 )  4 


COMMON/COMMOD/PARM.OTM.Q.ZOINIT.ZMODI  I  2 )  , ALF A  1  . BET A  1  , 
K  I  PATH! 20 )  , KPATH ( 20 )  , P I  . P2 , P 3 , PH . P5 . P6 . P7 . PB . 

K  P9.P10.P1 1 .P12.PI3.P1H 

C------------------- 

IF  < ZMOO ( 9) . GE . 1 . )  00  TO  202 


c-  -------  __________  PARAMETER  SETTING  .  10 

Pl-5.  n 

P2-0.6  13 

P3-0 . 80  13 

PH-  .  1 0  ]  4 

PS-2 . 5  15 

P6-0.9  16 

P7-1.5  17 

P8-.6  18 

P9-.3  19 

P10-.50  30 

PI  1  - . 7  2  1 

P 1 2-0 . 2  22 

P 13-22.  23 

PI  4- 1  . 3  24 

202  Continue  25 

ZP-ZM0D(Il*(T-ZM0D(6)l*ZM00(2l  26 

IF  ( ZM . GT . 0 . )  JO  TO  490  27 

C- . ALTITUDE  REPORT  MISSING:  COAST  TRACK  -  -  -  -  2B 

ZMOD ( 1 ) -ZP  29 

C  .  TO  BE  ADDEO:  IF  T-ZM00C31  GT  THRESHOLD.  DROP  TRACK  -  -  -  -  30 

GO  TO  802  31 

490  CONTINUE  32 

0ZM-ZM-ZM0D(41  33 

OBINS-ABS(OZM1 /Q  34 

ZMODI 9) -ZMOOI 9 1 »4 .  35 

IF  (ZM0DC9) .GT .P131  GO  TO  495  36 

C  .  -START-UP  SMOOTHING  --------  37 

ZMOD ( 1 1 -ZP+P9* ( ZM-ZP I  30 

ZM0D(2l-ZM00(2l*P4*(ZM-ZP)/iT-ZM0D(6))  39 

IF  (OZM.NE.O.)  00  TO  590  40 

GO  TO  802  41 

495  CONTINUE  42 

IF  (OZM.EQ.O.I  GO  TO  701  43 

C-  --------  -alt  report  changed  -  -  44 

ISGN-INTCS!GN( I . .DZMI I  45 

QSIGN-Q* ISON  46 

TEST-ZM0DI2) »DZM  47 

IF  ( TEST . LE .  1 00 . ANO . DB I  NS . EQ .  I  .  )  00  TO  580  48 

C . -----  -  -  TREND  EXISTS  - . -  49 

TPREV- t T-ZMOD (5) ) /OBINS  50 


Fig.  C.l.  FORTRAN  subroutine  used  to  update  a  track  at  4.7  second  intervals. 
This  routine  is  called  once  each  scan  following  the  scan  of  initialisation. 

As  a  programing  convenience,  the  setting  of  parameters  has  been  incorporated 
into  this  subroutine. 
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C  -  MO  TRANSITION  LOGIC  ----------  101 

701  CONTINUE  102 

ZMOD! 1 > *ZP*P9* IZM-ZP)  103 

TCUR*T-ZM00(5>  +  DT  10t 

TNOEX- ( TCUR-ZMODI 7 ) ) /OT  105 

IF  (ZMOO f 7 > . LT . 0 . 8*DT >  TNDE X - T CUR / ZMOO t 7 >  106 

IF  ( TNOEX . GT . P5 )  GO  TO  630  107 

IF  ( TNOEX . GE . P6 )  GO  TO  610  108 

c---. -------  NORMAL  UPDATE  -  NO  TRANSITION  EXPECTED  -----  109 

C------------  -rate  DECAYS  ACCORDING  TO  P3  .  110 

IF  ( ZMOO (SI. GE . 1 . I  GO  TO  60S  111 

ZMOD (2 > -ZMOD ( 2 ) *P3  112 

ZMOO ( 7 ) -0/  l ABS I ZMOO < 2 )  1  * .  1  )  113 

GO  TO  802  1  1  4 

610  CONTINUE  115 

C  -  -  -  -  BIN  OCCUPANCY  LONGER  THAN  EXPECTED  -  ALTER  EXTERNAL  RATE  -  -  -  116 

ZMOD IS) -SI GN(Q. ZMOO (2 ) ) / ( ZMOD (7 ) +(.H*ZM0D(7)+1.0«DT>*  117 

K  < TNDEX-0 . 41 • *2)  118 

ZMODI8 ) -AMAX 1  (2 .  . ZM00(8) -  1  .  1  119 

GO  TO  802  120 

630  CONTINUE  121 

C-- . TRANSITION  TO  LEVEL  FLIGHT  -  -  -  128 

ZMOD I  1  1 -ZM  1 2  3 

ZMOD (21*0.  124 

ZM00I7I-99.  125 

ZMOOIBl'O.  126 

ZMOO (101-0.  127 

i 32  CONTINUE  128 

IF  (ZM.GT.O.)  ZMOO ( 3 1  *  T  129 

ZM0D(6)-T  130 

RETURN  131 

END  132 

Fig.  C.l.  Continued. 
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0ELT-TPREV-ZM00I7)  51 

DZIO-OELT/OT  5 2 

IF  < ZMOD I 7 )  . UT . OT )  0 Z 1 0 • 1  . • l 0 T / ZMOD ( 7 I -OB  I  NS )  53 

IF  (ZHOO(S)  .LE.O.  )  GO  TO  540  5* 

IF  (ABSIOZIO)  . GT . P7  >  00  TO  540  55 

ZMODI 10)»P10»ZMODI 10I+DZ10  56 

IF  I  A BS I ZMOD I 10* ) . L£ . P 1 4  I  GO  TO  530  57 

c__.  ------------  EXCESS  RESIDUAL  DETECTED-  -------  58 

BETAI-Pll  59  • 

ZM0DI81-2.  60 

ZMOOI I0l«SIGNIP12, ZMOOI 10) )  61 

CO  TO  533  62 

530  CONTINUE  63 

C  -----  -  CURRENT  RATE  WITHIN  LIMITS  :  SMOOTH  -------  64 

Z7-ZM0DI7)  65 

BL l M* I Z7- l  .  )»*2/IZ7»*2*64.  )  66 

0ETA1 -AM AX l ( 1 . / (ZMOOIBl ♦ .6). SLIM. .10)  67 

ZM0018) -AMINl  IZM0DI8)*!  .  .  10.  )  68 

533  CONTINUE  69 

2MOO ( 7  I  - ( ZMOD ( 7 ) +BETA1 • ( OB  I  NS • TPRE V- ZMOD I  7 )  )  ) /  l 1  . +BET A  1 • I DB I  NS- I  .  )  70 

K  )  71 

ZMOD(2)*OSIGN/ZMOD(7)  72 

ZMOD I  1  ) -ZP  +  P9* I ZM-ZP I  73 

GO  TO  590  74 

540  CONTINUE  75 

C-- . . . REINITIALIZE  RATE  ------  76 

Z7PREV-ZM0D ( 7  I  77 

ZMOD l 7 ) • l . 2*TPREV+ . 05*0T  78 

ZM0D(2I-0SIGN/ZM0D<7)  79 

ZMODIIOl-O.  80 

ZMOD I  8  )  «  I .  81 

ZMOOI I ) -ZM-QSIGN/2+ZM0DI2) *DT/2  62 

GO  TO  590  83 

580  CONTINUE  84 

c. ----------  SET  TO  A  PRIORI  VALUE  FOR  SINGLE  TRANSITION  -  85 

ZMOOI 2 ) *P1 • I SGN  66 

ZMOOI  1  I *ZM-QS I GN/ 2* ZMOD 12 ) *DT/2  87 


ZMOD17)»QSIGN/ZMOOI2) 
ZMOD 181*0. 

ZMOD I  1 0 ) -0 . 

590  CONTINUE 


C_ -----------  UPOATE  LEVEL  TRANSITION  OATA  ------  92 

ZM00I4I-ZM  93 

ZMOD  I  5  I  •  T  9,4 

ZM0DI9) «ZMOD I  9) ♦ 1 0  .  95 

IF  I ZMOOI 3 ) . GE . ZM0DI6 ) )  00  TO  B02  96 

C  -  -  -  CORRECT  TRANSITION  TO  FALL  WITHIN  PERIOD  OF  MISSING  DATA  -  -  -  97 

ZMOD 15) «T*P8» I ZMOD I  3) -T»DT I  98 

GO  TO  802  99 


100 


Fig.  C.l.  Continued. 
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